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The free outflow hydrograph is studied in the case of a sudden water release from a 
reservoir. The outflow hydrograph is called free when it is not affected by the tailwater 
levels. Both the effect of the steep negative wave, created by sudden water release, and the 
effect of the flow resistance are analyzed through the use of a fictitious inflow hydrograph. 
The water accelerated by the steep negative wave movement, whether or not the flow resist- 
ance is taken into account, gives this fictitious hydrograph. The procedures are given for 
the computation of the wave celerities and the new water velocities along the reservoir, for 
the computation of wavefront heights and the maximum outflow discharges, and for the 
determination of the fictitious and total inflow hydrograph. The examples are given for the 
pyramidal reservoirs of rectangular and parabolic cross sections. A general procedure for 
the determination of the reservoir free outflow hydrograph in an approximate form is de- 
rived. 



1. Introduction 

1.1. Definition of the Free Outflow Hydrograph 

If the outflow hydrograph from a reservoir is un- 
affected by the tailwater levels or if such influence 
is negligible, it will be called the free outflow hydro- 
graph. In that case both the outflow hydrograph 
and the downstream propagation of the wave cre- 
ated by the sudden release of water can be analyzed 
separately, first the hydrograph and then the wave 
propagation. 

The reservoir free outflow hydrograph in the case 
of sudden water release depends on: (a) The area 
and shape of the outlet; (b) the type of flow (free 
surface flow or flow under pressure); (c) the shape 
characteristics of water body (storage function, cross 
sections, etc); (d) the inflow hydrograph; (e) the time 
interval in which a full opening of the outlet is 
created; (/) the movement of negative waves, espe- 
cially of steep negative waves; (g) the resistance to 
flow along the reservoir. 

The free outflow hydrograph can be approximated 
by taking into consideration only the most influen- 
tial of these factors, (a) to (d). Any classical reser- 
voir routing procedure will suffice for the computa- 
tion. The analytical solutions are possible in cases 
when the outflow rating curve, the storage function 
and the inflow hydrograph are given and can be 
approximated by simple expressions mathematically 
tractable [3]. 3 In that case the opening time is 
assumed as zero and the effects of steep negative 
wave movement and of flow resistance are neglected, 
assuming that at the very moment of fast opening 
the water in the reservoir was instantaneously ac- 
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celerated to give maximum discharge. This assump- 
tion is the source of the errors in the shape of the 
free hydrograph, if the steep negative wave is suffi- 
ciently high, as when the opening is large in compari- 
son with the reservoir cross sections. 

1.2. Purpose of the Study 

This study seeks the development of a simple 
procedure for the computation of free out (low hydro- 
graph, when the movement of the steep negative 
wave and the ilow resistance are taken into considera- 
tion. The procedure must be suitable for use with 
the available data, to the accuracy demanded of the 
resulting hydrographs, and it must require a reason- 
able computation time. 

1.3. Scope of the Study 

This study will investigate the influence of the 
steep negative wave and of flow resistance on the 
shape of the free outflow hydrograph. The knowl- 
edge of the shape of reservoir (cross sections) and of 
its roughness coefficients is needed. Since the open- 
ing time is very short in comparison with the outflow 
time of the major bulk of water, it can be taken as 
nearly zero. 

1.4. Basic Approach to the Problem 

The shape of the free outflow hydrograph, com- 
puted by neglecting the effect of the steep negative 
wave and of flow resistance, is to be considered as the 
first approximation of the true outflow hydrograph. 
This first shape is sufficiently accurate for small nega- 
tive wave heights in many practical cases, because 
the accuracy of background data, the routing pro- 
cedure and the accuracy demanded do not justify 
more precise computation. The first approximate 
shape of the free outflow hydrograph is thus con- 
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ceived as a gradual emptying of the body of water, 
starting from already accelerated water [3]. If the 
neglected factors were taken into consideration and 
if more data about the reservoir were available, the 
subsequent approximations would be closer to the 
true hydrograph. 

Two approaches can be used in the computation 
of the subsequent approximations of the true outflow 
hydrograph: (1) To develop a special procedure for 
computing the hydrograph directly from all given 
data and (2) to correct the previous approximate 
hydrograph by introducing the effects of new factors. 
In the second approach the effects are represented 
by a new fictitious inflow hydrograph. The varia- 
bles are time, the starting reservoir level, the charac- 
teristics of the steep negative wave, the length of 
reservoir, the flow resistance, etc. 

The new approximation of the true hydrograph 
can be computed from the first approximation by 
simple transformation using the fictitious function 
of inflow. Next, the sums of the real inflow hydro- 
graph and fictitious inflow hydrograph can be made. 
The new inflow hydrographs can be applied directly 
for the computation of the other approximations 
of the true hydrograph. This procedure will be used 
in this study. 

1.5. Effect of the Steep Negative Wave and of Flow 
Resistance on the Free Outflow Hydrograph 

The outflow from a reservoir in the case of sudden 
water release does create a steep negative wave which 
travels upstream along the reservoir. 

Every increase of outflow discharge also increases 
the head losses along the body of water. 

When a steep negative wave is created, its move- 
ment is influenced by: 

(a) The change of the wave area A x and of its 
shape (Ai is the vertical projection of the wavefront) 
and the change of the total cross section area A or 
of the area under the wave (A — Ai), figure 1, which 
change the propagation celerity; (b) the change of 
the direction and the width of reservoir, of the 
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Reservoir outflow hydrograph in case of sudden 
water release. 

(A) Schematic reservoir with negative wave movement; (B) two outflow 
hydrographs: 1— hydrograph with already accelerated water to Omax at the 
moment of release; 2— hydrograph greatly influenced by the steep negative 
wave. 




Figure 1. 



I. Steep negative wave 
2. Very gradually varied negative wave 

h= front height of steep wave 

C| ,C 2 ,C X = celerities of water body at rest 

H p H 2 ,H X = mean depths 

Figure 2. Two types of negative waves in case of sudden 
water release. 

Steep negative wave and long negative wave very gradually varied. 

shape and the roughness of its shores, which cause 
and affect the reflection and the refraction of the 
wave; (c) the resistance to the acceleration of water 
due to the movement of the wave, and the resistance 
of the accelerated water moving downstream between 
the wavefront and the outlet, which decrease the 
volume and the kinetic energy of the accelerated 
water. The flow resistance decreases the outflow 
discharges in the first part of the hydrograph and 
increases them in its tail part. 

As the propagation celerity of the negative wave 
gradually decreases in the upstream direction, 
because of the gradual decrease of mean depth in a 
reservoir with sloping bottom, the accelerated dis- 
charge per unit of time thus correspondingly de- 
creases. As the level at the outlet would be higher 
at any previous moment than the level corresponding 
to the accelerated arriving discharge at the next 
moment, a new level adaptation must be made to 
satisfy the continuity equation by a lowering of the 
level at the outlet, which causes the additional 
acceleration of water in the reservoir. The volume 
of water involved in this lowering is taken from 
surface layers of the body of water between the 
outlet and the steep negative wave. Therefore, a 
continuous propagation of small negative waves, 
changing very gradually, takes place from the outlet 
upstream. For each time unit a small drop in level 
at the outlet is created, wh ich propagates upstream 
with the celerity C=^gH x where H x is the mean 
depth (fig. 2). As EL, near the steep negative wave 
is lower than the mean depth of the lowest part of 
steep negative wave, the small waves can not reach 
the steep wave. As the process is continuous, the 
wave movement can be conceived as a group of very 
small waves making a new long wave, very gradually 
changing (fig. 2), with H x decreasing. 

1 .6. Approximations of the Free Outflow Hydrograph 

The combination of the movement of two types 
of negative waves, of a very steep wave (depression) 
and of a very long wave corresponding to gradually 
varied unsteady flow (with quite opposite char- 
acteristics), the complex reflections and refractions 
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of waves in a natural body of water, and the re- 
sistance to the flow make the exact computation 
of the free outflow hydrograph very complicated. 
Therefore, the solutions must be obtained in ap- 
proximations of a different degree of accuracy. The 
first approximation neglected the movement of 
negative waves and the flow resistance (as in [3]); 
the second approximation will include the addition 
of the steep wave effect; the third approximation 
will take account of the additional effect of flow 
resistance. These two approximations are treated 
in the subsequent part of this paper. Further 
approximations can be considered when practically 
feasible. 

2. Effect of the Steep Negative Wave on the 
Free Outflow Hydrograph 

2.1. General Assumptions 

As the breach of a dam or an accidental opening of 
some types of large gates are relatively fast processes 
(a matter of few seconds), the negative wave created 
by sudden water release is a steep wave traveling up- 
stream. Quite different is the case when the opening 
of outlets is very slow, the outflow rating curve 
changing with time as the opening becomes larger. 
Q is a function of opening dimensions and of time. 
This type of opening will not be studied here. 

The steep negative wave undergoes three types of 
changes as it travels upstream: (a) Reflections and 
refractions; (b) decrease of propagation celerity due 
to the change of cross section; and (c) loss of energy 
due to the resistance. 

The reflection and refraction of the steep negative 
wave are complex because of the factors which es- 
sentially change continuously along the reservoir. 
The process is, however, somewhat less complex in 
channels of nearly regular shape. The effects of the 
reflection and refraction of the steep negative wave 
will be assumed in (he second and third approxima- 
tions as negligible for the valley reservoirs. 

The effect of the negative wave ends when its prop- 
agation celerity C in water at rest equals the menu 
flow velocity V already existent (C=V), or when the 
wave arrives at the shallow regions, if water is in rest 
during the wave movement. 

The effect of resistance will be neglected in the sec- 
ond approximation, which is justified especially for 
large bodies of water and moderate heights of steep 
waves. It will be considered in the third approxima- 
tion (sec. 3 of this paper). 

The main factor in this second approximation is 
the additional influence of the existence and move- 
ment of steep negative wave. 

The flow through reservoirs or lakes, where the 
mean velocities are relatively low and the cross sec- 
tions change gradually, is a gradually varied flow, 
either steady or unsteady. 

Before the analysis of the effect of the steep nega- 
tive wave on the outflow hydrograph is begun, two 
problems will be studied: (1) The propagation of the 
wave along the water body and (2) the steep wave 
height and the outflow discharge shortly after the 
sudden water release. 




Figure 3. Schematic representation of the transformation on 
the front shape of steep negative wave, due to different celerities 
of wave elements. 



The shape of a steep negative wave moving along 
the body of water is not stable (fig. 3). The reason 
is that the parts of wavefront travel with different 
celerities. The steepness of the wave decreases with 
the distance traveled, and the final shape will be ob- 
tained when the slope of the wavefront is essentially 
equal to the slope of the bottom. Because of these 
differences in celerities of different points of the wave, 
it is necessary always to define C (as that of the high- 
est point, of the lowest point, of the center of gravity, 
of the mean height of wave, etc.). The approxima- 
tions in this study deal with the schematic negative 
wave, considered as a discontinuous steep wave. 

2.2. Mean Celerity of the Steep Negative Wave in 
a Cross Section 

The steep negative wave arrives at a cross section 
of the body of water, X distance from the outlet, 
at time t from the opening time, with wavefront 
height h. The celerity of the wave 4 in water at rest 
is C at X. The mean velocity of water because of 
inflow is V. The cross section characteristics are 
(fig. 1): Area A before the wave arrival; area A x 
of the wavefront height; distance/ of the center 
of gravity (G of area A x from the new level (or from 
the starting level f x = h—f); area A—A h after the 
wave has passed. The effect of the bottom slope 
is indirectly introduced by the change of A, A ly 
jj u and^l-A. 

Two equations, continuity and momentum, give 
the celerity of the steep negative wave and the 
increase of the mean water velocity in a cross 
section of reservoir, due to the acceleration caused 
by the wave movement, as functions of section 
characteristics [1, 2]. For a small time interval dt 
the celerity can be considered constant in the vi- 
cinity of the cross section at X, regardless of (he 
bottom slope which has little influence on the cross 
section characteristics for short path length (fig. 4). 

The steep negative wave travels upstream a 
distance (C-V)dt during the time dt (fig. 4), and 
the water travels downstream the distance (V-\- U)dt. 
An observer, traveling at the section X with 
the velocity V downstream, after time dt, will be 
at the section X— Vdt, or a distance Vdt downstream. 
To him the distance traveled by the wave in the 
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Volume 1-2-3-4-1 = A,C x At 
= Volume 4-5-6-7-4 = (A-A,)U X At 



A, C X =(A-A,)U X 



Figure 4. Scheme for the derivation of continuity equation in 
case of movement of steep negative waves. 



upstream direction is Cdt and the distance traveled 
by the water downstream is Udt. 

The continuity equation gives for the observer 
(or also in case V=0) 

A 1 Cdt=(A-A 1 )Udt, 

so that generally for any X 



c 



=<-> 



a) 



where C is the propagation celerity of the steep 
negative wave in water at rest and C— V is the prop- 
agation celerity in moving water with velocity 
V (V<^C). Equation (1) can be obtained also by 
equating the total water discharge entering the 
small volume Adx, to the discharge leaving that 
volume 

VA+(C-V)A 1 =(V+U) (A- A,). 

For the momentum equation, it is assumed that 
the acceleration of water from the velocity V to the 
velocity V-\-U of the volume of water (A—A^Cdt, 
in time dt and at the section X (figs. 4 and 5) is 
made without loss of energy. 

The volume element (A— Ay) Cdt, which is accele- 
rated from V to the velocity V+ U between times 
t and t-\-dt, has a momentum (A—A^CUdt and a 
rate of change of momentum (A—A^CU. The 
force F which creates this change is the difference 
of pressure on the areas A and A—A u (fig. 5), 
F=gh(A—Ai)+gfA, with/ already defined. 

This force is equal to the rate of change of momen- 
tum, so that 



CU=gh+gf 



A { 



A-A, 



The celerity C at any section in water at rest, 
after using eq (1), is 




A= area l-2-345i 
A, = area I -2-3-5-1 
A-A,«area 3-45-3 

Figure 5. Scheme for the derivation of momentum equation in 
case of steep negative wave movements. 



where a=A/A l . The ratio a>l, and theoretically 
can be between the limits 1 and infinity. If h 
increases, (a— I) decreases faster than h and /, 
so that the mean celerity of the wavefront decreases 
with the increase of h. With the mean depths of 
areas A and A ly H m =A[B and h m =A 1 /B (fig. 5), 
with B the width at the surface level, and putting 
fi = h—f, eq (2) becomes 



°=aK*-kt''> 



(3) 



When h/hm is nearly 1 and f x can be neglected 
in comparison with H m (for small values of h) then 



C=ygH m , 



(4) 



which is the known expression for the propagation 
of small surface disturbances in shallow water, with 
the mean depth H m . 

The increase of the mean water velocity is, there- 
fore, from eqs (1) and (2) 



U=— x jgh(a-l)+gf, 



(5) 



or from eqs (1) and (3) 



C=jgh(a-l)+gf, 



(2) 



u =e^kM h 4-i)> ^ 



which approaches zero when h approaches zero. 

The discharge P W =U(A—A 1 ) is the accelerated 
discharge by the movement of the steep negative 
wave, which will be considered as the fictitious inflow 
in the reservoir. 

As H m , h m , and j x are functions of cross section 
characteristics and of height of wavefront h, both 
C and U are functions of those parameters. The 
cross section characteristics are functions of the water 
level. For a given section, level, and height h the 
values C and U are uniquely determined. 
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2.3. Computation of the Initial Front Height of Steep 
Negative Wave and of the Maximum Outflow 
Discharge 

If an opening in a reservoir is created with constant 
area and shape, its hydraulic characteristics are 
defined by the outflow rating curve, generally of the 
type Q=bk T , where b and r depend on the area, shape 
and position of the opening, as well as on the type 
of flow, with k the total energy head above the level 
of the outflow zero (fig. 1). The opening height, as 
he difference between the level at the moment of 
opening and the level of outflow zero, is h . It is 
the sum of two heights: h, the front height of the 
steep negative wave and k , the difference of the 
level just upstream from the outlet (where the 
depression due to the outflow is negligible) and the 
level of outflow zero, so that h =h+k . The energy 
head is 

a(v±uy (7) 



k=k 



2g 



where V is the velocity near the opening before it is 
created and a is the kinetic energy correction factor. 

For the small time interval, At, after the opening, 
the water flowing through the opening is the sum of 
the inflow discharge P=AV and the accelerated 
discharge due to the moving steep negative wave. 
The interval At must be sufficiently long for the 
negative wave to develop but not so long as to 
increase the outflow by the further significant drop 
in water level of the stretch of reservoir between the 
outlet and the steep negative wavefront. 

The continuity equation between supply and out- 
flow gives 



Q =AV+A 1 (C-V)=bI^ t 



(8) 



which can be written as [4, p. 213, without last term] 

Q =AV+A 1 (C- V) =b (ko+ a(V ^ g Ur J- (9) 

For small value of velocity head in comparison 
with k , eq (9) can be expressed in the approximate 
form 

Q =b^=bkl[lV^^] (9a) 

The second term of [4, p. 213] in form of (T^/2#) 3 / 2 , 
which takes here the form 



4^FI 



(9b) 



and should be added on the right side of eq (9), can 
be neglected for reservoir outflows in comparison 
with the first term. The celerity C is given by eq 
(2) or (3), and velocity U by eq (5) or (6). 

For given h , b, and r (values or functions depend- 
ent only on the outlet characteristics), and given 
V and a, the magnitudes A u C, U, and k =h —h 



are functions of h and of the cross section character- 
istics near the outflow, at distance AX=(C— V)At. 
When the cross section at distance AX from the 
outlet is selected, eq (9) enables the computation 
of the value h. As the discharges around the maxi- 
mum discharge change slowly (fig. 1), the selection 
of AX is easier. It is, therefore, sufficient to select 
an average cross section some hundred meters 
upstream from the outlet. For given h, the values 
k and U, or A x and C can be obtained, so that eq 
(9) also gives Q , the discharge outflow at the time 
At (which is the maximum outflow if AX is well 
selected). 

The computation of h is to be made by successive 
approximations until eq (9) is satisfied for the com- 
puted h. Numerical or graphical procedures are 
feasible for that computation. 

A simple graphical procedure for the determination 
of h, Qq and downstream tailwater level H d is shown 
in figure 6. Three curves as the functions of h and 
independent of the height h of the opening are first 
plotted: 

Q 1 =A 1 (C-V)+AV, 



Ah 



and 



= a(V+U) 2 

Q*=KH d ), 



where Q 3 is the rating curve of the tailwater channel. 
For the given opening height h the fourth curve 
Q 2 =bk T is plotted. For some types of outlets, the 
coefficient b is a function of k /(H —h ) and of some 
other parameters, so that approximate values must 
be first introduced. For convenience the velocity 
head Ak is plotted with the same scale as h , h, k, 
and k . For an assumed value h or k (point 1) the 
discharge Qi = A x (C— U)-\-AV is obtained (curve 1, 
point 3), as well as Ak (curve 3, point 2). By adding 
Ak to k (point 4) the value k is obtained, and 
Q 2 =bk T (curve 2, point 5). The adjustments have 
to be made for h until Qi = Q 2 > If ft = $2, the value 
of h is obtained. From that value the discharge 



(J) Q,=A,(C-V) + AV 
(D Q 2 ~- bk r 




k=K,4 



alV-HJr 



® Oy t (H d ) 



H m =mean depth upstream 
of the outlet 



Figure 6. Graphic procedure for the determination of front 
height of steep negative wave, of maximum discharge, and of 
maximum tailwater level. 

1— accelerated discharge plus inflow Po as function of wavefront height; 2— 
outflow rating curve; 3— total velocity head as the [unction of wavefronl height; 

4— tailwater rating curve. 
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Qo is determined (point 6). For a given Q (curve 4, 
point 7) the maximum downstream level H d (point 8) 
is obtained. When H d </3(H —h ), the outflow is 
uninfluenced by tailwater level. The coefficient 
&> 1 is a correction factor which takes into account 
that, due to the effect of kinetic energy of outflow 
water, the tailwater levels somewhat higher than 
(H —h ) do not influence the outflow. 

2.4. Relation of the Height of Steep Negative Waves 
to the Opening Height of a Rectangular Outlet 

For the given opening, the values h and Q depend 
on the shape of three curves, 1,2, and 3 of figure 6. 
The accuracy of curve 2 (Q=bk r ) depends on the 
accuracy of the dimensionless outflow coefficient m 
which enters into the expression for the coefficient b. 

The case of a reservoir, with water at rest, of 
approximately rectangular cross section A=B H m 
near the outlet will be treated here. The bottom 
of small length of reservoir AX can be assumed 
horizontal, especially under natural conditions with 
deposited sediments. The simple case of a rectangu- 
lar sharp-edged opening will be analyzed here as an 
example, with area Bh , where h is much smaller 
than H m (fig. 6), and B<B . Equation (9) will be 
used, and according to [4] and [5], with different 
symbols used here, 



can be written as 



m= 



Q 



Q 5 _ b B 



or 



Bhjgh B^gh 3 ' 2 B B^ g B 

— Hah,'*'™) <10) 



with k =h —h and R and W the Reynolds and 
Weber numbers, respectively. 

In this expression B is taken in the denominator 
instead of the usual B, so that the right side of eq. 
(10) is multiplied by B /B, apart from the contraction 
effect dependent on B/B , in order to facilitate the 
solution of eq (9), regardless of the changing width 
B of the opening. 

Putting fe /£L=s and h/H m =p, the ratio h/h =p/s, 
of the wave height h and the opening height h will 
be related to s. For the rectangular sharp-edged 
opening 



m- 



Bo 
' B 



f(&'T3*> w } (11) 



where m is B /B times higher than the usual coefficient 
m. Other functions of m for other shapes of outlet 
have to be used, depending on the outlet. 

The coefficient m can be thus expressed as a func- 
tion of dimensionless numbers, so that for each outlet 
and its characteristics the special relation of eq (11) 
or of other type can be obtained. For given B, B , 
H m and h 0j the coefficient m can be expressed as a 
function of h, or of p = h/H m . 

Equation (9) under the above assumptions and in 
terms of the dimensionless numbers introduced here 



s=p+ 



r) 2 /3 



/// 



2/3 



o-r 



a p 2 (2-p) 
4 (l-2>) 2 ' 



(12) 



s are more 



If the function of eq (11) is introduced in eq (12), 
the relationship p=f(s) is obtained. For the practi- 
cal solution it is suitable to take m as a parameter 
and to develop a family of curves for eq (12). Figure 
7 gives this family of p=f(s) for m=0.25, 0.50, 0.75, 
1.00, 1.25, and 1.50, and figure 8 gives p/s=f(s) for 
w=0.25, 0.50, 0.75, and 1.00. The families of curves 

m== /(^) or m>=f ( ) for given values 

convenient to use. Both are given in figures 7 and 8 
as dashed lines. For these curves a=1.05 was taken 
as the most applicable constant value. 

When eq (11) is determined for a particular case 
of rectangular opening, with s=h /H m known, the 
intersection of the curve of eq (11) and the curve 
m=f(h/H m )=f(p) or rn=f(h/h )=f(p/s) from figure 7 
or 8 for that s, gives the solution either for p=h/H m 
or p/s=h/h . 

If the velocity V in the body of water near the 
outlet, before the sudden release of additional water, 
can not be neglected in comparison with the new 
velocity (F+?7), it must be introduced in dimen- 
sionl ess form in eq. (12). In that case, for w= 
V /^gH m = V / C my with C m =^[gH~ m , and V Q =P/B H m , 




Figure 7. Relation of the relative height of wavefront (p=h/H m ) 
to relative opening height (s = holH m ) where H m is mean depth 
of reservoir near outlet, with the outflow coefficient m as 
parameter, for rectangular opening; the second family is the 
relation of m and p with s as a parameter. 
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Figure 8. Relation of the ratio of wavefront height and opening 
height (h/h ) to relative opening height (ho/H m ) with outflow 
coefficient m as parameter for rectangular opening; dashed 
lines represent the relation of m and p/s with s as parameter. 



(13) 



For each given value of w a family of curves 
m=f(h/H m ) or 7n=j(h/h ) can be developed. 

Equations (12) and (13) in parametric form can 
be developed for any regular-shaped outlets (rec- 
tangular, triangular, semicircular, circular, trape- 
zoidal, parabolic, semielHptic, etc.). With s, or with 
s and w, as parameters, the families of curves 
7 n=f(h>/K) or m=f(h/H m ) can be developed once 
for all cases. With these families available, the 
determined real function m=J(h/h ) for a particular 
outlet of given regular shape will be sufficient for 
the determination of the height h of the steep nega- 
tive wave, and of the maximum discharge of outflow 
hydrograph Qo by using eq (8). The nomographs 
can also be developed for the relations in eqs (12) 
and (13). 

If m is selected to be only the function of h, then 
eq (12) would give a stable relationship between p 
and 6 j , or eq (13) would give a family of curves with 
w as parameter. 

2.5. Procedure for the Computation of a Free 
Outflow Hydrograph 

In the case of no inflow into the reservoir, the 
celerity C depends on the wave height h and on the 
depth II of reservoir at the outlet before the sudden 
water release. The celerity is, therefore, the general 



0=f ( x >ww) 



(14) 



where x=l — X/L is the relative distance of a cross 
section from the upstream end of reservoir, X is 
the distance from the outlet, L is the length of 
reservoir, H is any possible reservoir depth at the 
outlet, H is the maximum operational depth, h 
the height of the steep negative wave. 

For a given relative parameter H/H , a family of 
curves C=f(x, h/H ), can be plotted. 

For a large reservoir, H and h will be assumed as 
the initial or boundary conditions. For a given 
cross section x the values A, A u a, /, C, and U, 
entering into eqs (2), (3), (5), and (6), are also 
given. 

The general procedure for the second approxima- 
tion of taking into consideration the effect of the 
movement of steep negative wave in the case of 
valley reservoirs is: 

(1) To compute the wave height h for given water 
depth //at the outlet, and to determine the maximum 
outflow discharge Q , as is shown in subsections 2.3 
and 2.4, eqs (7) to (13), figures 6 to 8. 

(2) To verify that the outflow hydrograph is free, 
or that the tail water level is lower than the maximum 
tailwater level which does not affect the outflow 
(fig- 6); 

(3) To compute the celerity C along the body of 
water out of given cross section characteristics (A, 
A { and/), eqs (2) and (3); 

(4) To compute the velocity V along the body of 
water for the inflow / > dining the wave movement, 
it there is an inflow; 

(5) To determine the difference (C—V) along the 
body of water; 

(6) To compute the accelerated discharges by 
steep negative wave movement along the body of 
water P w i = CA 1 with inflow zero, or P w2 — {C—V)A X 
with inflow T P; 

(7) To determine the time, when the wave arrives 
at the section X and when the acceleration takes 



C x dx 

place there, by £= I -^ — ^f 
Jx=o k — V 

(L= length of reservoir) by 

, T C l dx 



or 



for 1-f- 



(15) 



with £ the total time of wave travel from the outlet 
to the upstream end of reservoir (from x=l to 
3=0); 

(8) To find the relation between the accelerated 
discharge and the time of acceleration, P w =J(t). 
For given x the pair of values P w and t gives P w =j{t). 

The integration of P w over time from £=0 to 
£=£o gives the volume S h of the reservoir layer h, 
determined simply from the storage function S=f(H), 
figure 10, which can serve to check the function 

*W(«). 
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If the inflow during time t cannot be taken as 
constant (i.e., the natural inflow hydrograph or the 
release of water from upstream reservoirs are very- 
unsteady), the velocity V at a cross section x is also 
dependent upon the time as are C— V, Ai and P w . 
The inflow can be either a steep or a gradual wave and 
the determination of P w =j(t) becomes more complex. 

If the inflow is in form of a steep positive or nega- 
tive wave, the procedure is similar to the above given 
movement of steep negative wave until two waves 
traveling in opposite directions meet. The resulting 
waves will then move in both directions. 

In the case where the inflow hydrograph P=j(t) is 
gradual^ varied, and where the steep negative wave 
moves faster upstream than the other negative waves 
created downstream of the steep wave, a resulting 
new inflow hydrograph P n =f(t), composed of P=j(t) 
and the accelerated discharges by steep negative 
wave movement, can be determined either by exact 
or approximate methods. The position of the wave 
P=f(t) along the reservoir has to be determined. 
Therefore, the velocity V for given places and the 
times at which the steep negative wave arrives at 
those places must be computed in order to obtain the 

celerity C— V. The celerity is now C=f ( x, t,jjj jj- ), 

with a new variable, t. 

The time t of wave travel from one end of the 
reservoir to the other is the factor which will deter- 
mine whether or not the inflow hydrograph can be 
replaced by a constant inflow P . Since t for the 
majority of valley reservoirs is measured in hours, 
constant inflow can often be assumed during the 
movement of steep negative wave. 

Due to the neglected factors — -resistance to the 
flow, reflection and refraction of the wavefront and 
transformation of the wavefront due to different 
celerities of its parts — -the computed function P w =f(t) 
is an approximation of the true relationship. 

The procedure for the second approximation of the 
free outflow hydrograph would involve the following 
values and functions : 

(1) The outflow head k=h+ a ^ V ^ ^ (2) the 

outflow rating curve Q=b k r ; (3) the characteristics 
of cross sections along the reservoir, and (4) the total 
inflow function P n =f(t), which gives for each time t 
the discharge P n =P+P w , and for the corresponding 
x, also, the place of inflow P n . The use of the two 
basic de Saint Venant equations for gradually varied 
unsteady flow will enable the computation of the free 
outflow hydrograph in cases where the continuity 
equation alone can not give reliable results. 

In valley reservoirs, where the height of the nega- 
tive wave is not too large, the use of the continuity 
equation alone, instead of both continuity and mo- 
mentum, gives an approximate free outflow hydro- 
graph, sufficiently accurate in the light of background 
data, assumptions and the needed accuracy of results. 

The higher the steep negative wave and the inflow 
P, the greater is the velocity (V+ ■U) l so that the use 
of the continuity equation alone will give less accurate 
results. Since the effect of neglected factors in steep 



negative wave movement is the highest in shallow 
waters (upstream part of reservoirs) and lowest in 
the deep waters (lower part of reservoirs) and since 
P w is low in the first and high in the second region, 
the use of only the continuity equation is generally 
justified. In that case the integration of the storage 
differential equation (continuity equation) of the 
type P— Q=dS/dt, with P replaced by the new inflow 
hydrograph P n =P-\-P w , can be made by a numerical, 
graphical, or any other appropriate procedure to ob- 
tain the free outflow hydrograph Q=f(t). 

Two simple cases will be taken as examples : (a) A 
pyramidal reservoir with rectangular sections (B x , 
H x ), with the bottom slope H /L, where H is its 
depth at outlet and L reservoir length, and with B 
the width of reservoir at the outlet; (b) a parabolic 
section of the type B=Bq tJH/H , where B and H are 
width and depth respectively of reservoir at the out- 
let, B and H any width and depth, with the bottom 
slope Hq/L, and B x and H x the surface width and 
water depth at a cross section X distance from outlet. 

2.6. Rectangular Cross Sections 

Figure 9 gives the simple shape of the reservoir and 
the resulting functions along it for the inflow P=0, 
The steep negative wave has the front height h with 
p=h/H , anda?=l— (X/L), so that B x /B =H X /H =x. 

The total area of any cross section is A=x 2 B H 
and the wavefront area Ai=xhB ; j=h/2, so that 
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Figure 9. Celerity, velocity of accelerated water, and accelerated 
discharges along a pyramidal schematic reservoir with rec- 
tangular cross sections and rectangular opening of the same 
width as the reservoir (upper figure) . 
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from eq. (2) 



C=C, 



'°\r _ 2' 



(16) 



with Co=^gH , the maximum celerity of the upper- 
most part of steep negative wave near the outlet. 
In the case of constant inflow P and V=P /A, 



C-V=G 






(17) 



with Vo=P /B H Q the velocity at the outlet before 
the sudden water release. The new velocity along 
the body of water after the wave passes is, according 
to eq (5), 



v+u=^+a 



x V 



(18) 



For small values of x eq (16), (17), and (18) have no 
practical meaning, because G and U are very high, 
and the resistance plays the major role, so that they 
should be used for values slightly higher than p. 
The accelerated discharge plus the inflow P is now 
P n =P Q +P , with 

P^hBofcoxJx^^y (19) 

The travel time from the outlet to a cross section 
X distance, in the case V=0, is 



. C x dX r ndx 2L n . 



(20) 



If x=0, the travel time along the reservoir for Vo=0 is 

2Z 



*0 = 



a 



(21) 



t 



and x=l— -=1 — T, where T=t/t . With tlie in- 
to 
crease of P the time U increases and the discharge 
P w decreases, though P n =P +P w increases. 

The integration of the £— function for Vo^O in 
the form 



< 



x 2 dx 



CqX' 



V a 



"2 Vo 



(22) 



is best made numerically. 

The functions C, U, andF w are given in logarithmic 
scales in figure 9 for the following values: H = 100 m; 
B = 300 m;h=10 m;C and U in m/sec and Q in m 3 /sec. 

2.7. Parabolic Cross Section 

For given h and H , and two cases of the inflow: 
P=0 and P=Pq, figure 10 gives the following func- 
tions depending on x=l — (XjL) for the parabolic 
cross sections of the type B=d^H H, where d=B /H ; 



A; A X ) a=A/A i; C (or C-V); U (or U+V); V=P /A; 
t= 77 ( or ^ = Ti — ~tt)> I j i ( or P2), accelerated 

discharge P w or P n =P +A l (C—V); and finally P x 
(or P 2 ) as a function of time t\ (or t 2 ). 

2.8. Shape of the Hydrograph in the First Moments 
After Release 

A few moments after the release, at the time in- 
terval At, a full steep negative wave of height h is 
developed (fig. 1). At that moment the outflow is 
maximum. The approximate solution for the shape 
of the rising part of the hydrograph after the sudden 
water release can be obtained as a straight line for 
an ideal fluid by using the two-dimensional hydro- 
dynamic equations of motion expressed in Lagran- 
gian form [6]. 

If the straight line rise is assumed for the total 
rise, the procedure for the determination of the slope 
of the hydrograph is simple. The total outflow of 
water from the time of opening to the time A£ is ap- 
proximately the volume of water confined between 
the outflow section and the section at AX= (C— V)At, 
where (C—V) is the wave celerity at the beginning, 
and between the starting depth of water H at the 
outlet and the new depth H—h. This volume is AS. 
If it is taken that the maximum discharge Q takes 
place at time A^<A^, instead of at A£, the continuity 
equation gives K(A£+A£ 1 )$ =AS', and A^ can be ob- 
tained (fig. 1). The mean slope of the straight line 
for the rising part of the hydrograph is then Qo/Ati. 
For small h, the values AS and AT are very small, 
and the rise of hydrograph can be taken as vertical. 

Figure 11 gives an example of the difference be- 
tween the outflow hydrograph of the first approxi- 
mation, curve 1 (emptying of the reservoir with the 
negative wave movement neglected) and the outflow 
hydrograph of the second approximation, curve 3 
(with the steep negative wave taken into considera- 
tion, but with the flow resistance neglected) for the 
Savage River Reservoir near Cumberland, Md. A 
rectangular opening (250 m by 18.5 m) with sudden 
water release in the case of a full reservoir is 
assumed. Curve 2 gives the fictitious inflow 
hydrograph computed by the method given above. 
Though the total volume of the two hydrographs, 1 
and 3, is the same, the peak discharge in the second 
approximation is smaller for about one-third. 

3. Effect of Flow Resistance on the Free 
Outflow Hydrograph 

3.1. Two Aspects of the Flow Resistance 

The effect of the flow resistance on the free out- 
flow hydrograph depends on the new velocities along 
the reservoir. The greater the wavefront height h, 
the less the remaining flow area A—A x after the 
steep negative wave has passed, and the greater the 
inflow and the roughness coefficient in the reservoir, 
the greater will be the flow resistance effect. 
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Figure 10. Schematic representation of the cross-section characteristics and h 
nitudes along a pyramidal reservoir with parabolic cross sections. 

(A) reservoir with areas A and Ai and their ratio a=A/Ai along the reservoir with a cross section at left and storage 
function of the reservoir at right; (B) celerity C and velocity U of accelerated water, due to the steep wave movement 
along the reservoir in case where water is at rest, and celerity C- V and velocity V=Po/A and *V+ U along the reservoir 
in case where inflow of water in reservoir is Pa; (C) accelerated discharge and times at which the acceleration take place 
along the reservoir for two cases: inflow is O and inflow is Po; (1) fictitious and total inflow hydro-graphs for reservoir in two 
cases: inflow is O and inflow is Po (third-dashed-curve gives fictitious hydrograph in case inflow P=0). 



Two aspects of flow resistance can be taken into 
consideration: (1) Head losses along the reservoir, 
which can be taken into account by introducing the 
friction slope into the second or dynamic de Saint 
Venant equation for unsteady flow gradually varied; 
and (2) decreasing of front height of the steep nega- 
tive wave, because of the downstream increase of 
level due to the resistance head as the wave pro- 
gresses upstream. Hence, the effective potential 
energy of wave (energy to be transformed into kinetic 
energy of accelerated water) will be smaller than the 
total energy corresponding to the initial height h. 

The first aspect of resistance can be taken into 
account when the routing procedure is based on the 
two de Saint Venant equations. It is, however, 
possible to take into consideration the resistance 
effect on steep negative wave movement when either 
the continuity equation of the type P—Q=dS/dt 
or both the continuity and dynamic equations are 
used for the routing. The effect of flow resistance 
on the fictitious inflow hydrograph, caused by the 
steep negative wave movement, will be dealt with 
in this investigation. 



3.2. Effect of the Flow Resistance on the Steep_Nega- 
tive Wave 

If a cross section at distance X from the outlet has 
the following values determined: A, A 1 , f, C, V, U, 
and if the resistance effect is introduced, the height 
h will decrease along the reservoir, with the slope 



dh_n 2 (V+U) 2 



dx 



R^ 



(23) 



where the Manning formula for the flow resistance is 
used in metric system, where n is the roughness 
coefficient and E the hydraulic radius of area 
{A—Ax). At the distance x=l — (X/L) the new 
wave height hi is now 



hi=h — h : 



-»-/, 



■■nW+UYdx 

#4/3 



(24) 



where h is the height near the outlet, hi the height of 
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Figure 11. Three hydrographs of the Savage River Reservoir for the rectangular opening 250 
m by 18.5 m: 1 — first approximation hydrography when the movement of alt waves is 
neglected; 2 — fictitious inflow hydrography when the water accelerated by the movement of 
steep negative wave is considered as the inflow; 3 — second approximation hydrograph, when 
the movement of steep negative wave is taken into account. 



wave at section x, and h x the total resistance head 
from the outlet to the section x. Since U and R are 
the functions of hi, the solution has to be made by 
successive approximations. As h x changes slowly, 
it is often sufficient to start with A=constant for the 
computation of U and R, and then correct values 
so obtained by successive approximations. 

In fact, the potential energy of the steep negative 
wave is not all transformed into kinetic energy of 
accelerated water, one part of it being left as potential 
energy to be used as the flow resistance head. 

The integration of the last term of eq (24) along 
the reservoir from the outlet to the place where 
C— V=0 would give the total resistance head. The 
volume S r which is not accelerated by the wave be- 
cause of the resistance is part of the total volume 
S h (fig. 10) corresponding to wavefront height h, and 



is 



S r = I BJijdx, 



(25) 



where a?o is the place where C— V=0 and B x is the 
mean width of the area element for total resistance 
head h x at the section x. 

The analytical integration of eqs (24) and (25) is 
not tractable even in simple cases like the pyramidal 
reservoir of rectangular cross section, figure 9, be- 
cause the values U and R for a section x depend on 
the sum of previous dh, and because of the expression 
Rv\ 



The practical procedure would be to compute U 
a ikI R based on constant value h along the reservoir, 
and to compute the slope of eq (23). By integrating 



h 



r x dh 

Ji dx 



dx the function of total resistance head 



h r along the reservoir is obtained. Then hi=h—h x 
gives the effective head of wavefront. With these 
values of It, new values [/and B are again computed, 
and eqs (23) and (24) can he obtained as a second 
approximation. Succeeding approximations can be 
made by continuing the same procedure. 

Due to the fact that the total velocity (V+ U 
increases in the upstream direction while the acceler- 
ated discharge decreases, it is usually sufficient to 
compute h x to the second or third approximation. 

With the new values h l = h—h x =j{x) it is now 
possible to determine a more exact fictitious hydro- 
graph P w . 

For the computed function h x =f(x) the width B x 
has to be determined since B x =f(x). The integra- 
tion of eq (25) gives S T . The product of BJh x along 
the reservoir makes possible the approximate compu- 
tation of the mean distance j r of total volume S r 
from the level H () —h: 

fr= Wr = ^r£ BxhldX ' (26) 

where E T is the potential energy corresponding to the 
volume S r and the resistance head and y is the unit 
weight of water. 
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The energy of the steep negative wave (fig. 10) is 



4. Discussion 



E w =yf Q S h 



(27) 



where S h is total volume of reservoir between depths 
H —h and H , and/ is the distance from the center of 
gravity of S h to the H —h (fig. 10). 
The ratio of energies 



E T = frSr 



(28) 



can be taken as a measure of the flow resistance 
effect on the transformation of the front height of 
steep negative wave. 

As h x increases somewhat faster than B x decreases 
in the upstream direction, B x h x increases from the 
outlet upstream. The effect of flow resistance 
decreases the accelerated flow P w , slower at its highest 
values near the outlet and faster at its lowest values 
at the upstream end of the reservoir. The fictitious 
inflow hydrograph will be now closer to a straight 
line. As a practical approximation for rapid compu- 
tations this hydrograph can be taken as triangular- 
shaped for a valley reservoir. 

3.3. Simplest Procedure for the Computation of the 
Steep Negative Wave Effect With Resistance 
Influence 

A simple procedure for the computation of the 
approximate free outflow hydrograph in this case 
would be: 

(1) Compute h and Q Q ; (2) compute A, A u j, C, 
V, U along the reservoir for constant h; (3) determine 
h x along the reservoir; (4) recompute all values under 
(2) ; (5) recompute h x and all values under (2) until a 
sufficient approximation of the true values is ob- 
tained; (6) compute B x as/(#); (7) compute S r for 
given h x and B x curves; (8) determine S h from the 
storage function; (9) determine the difference S h —S r ; 
(10) determine the triangular shape of fictitious 
inflow hydrograph P w with the maximum discharge 
Oo and the total time of fictitious hydrograph. 



U~- 



2(S h -S r ) 



Qo 



(29) 



The sum of the original and of the fictitious inflow 
hydrograph gives the total inflow hydrograph. The 
new inflow hydrograph, the maximum discharge Q 
at the new starting level k (fig. 6), the outflow rating 
curve and the storage function are now sufficient 
elements for any suitable routing procedure, based on 
the continuity equation only, of obtaining the free 
outflow hydrograph. 

If the routing procedure involves the use of both 
the continuity and dynamic equations for unsteady 
flow which is gradually varied, the computed fictitious 
inflow discharges P w , with their positions in time and 
along the reservoir, and the corresponding reservoir 
depths Ho — hi for the same times and places are the 
boundary conditions to be introduced in those 
equations. 



The computed outflow hydrographs from valley 
reservoirs in the case of sudden water release can be 
only approximations of the true hydrographs. The 
shape and the resistance to flow of valley reservoirs 
as well as the water movement itself along them are 
usually very complex. The closer an approximated 
hydrograph is to the true hydrograph, the more work 
is involved in the computation. The digital com- 
puters will, in some cases, give the most accurate 
results. The use of computers, however, is justified 
only in cases where very reliable and sufficient records 
and data about the reservoir, outlet, and inflow are 
available and where the great precision of results and 
speed of computation make their use economically 
feasible. Therefore, the use of simple numerical or 
graphical procedures for the computation of approxi- 
mate outflow hydrographs is often very convenient, 
sensible, and economically desirable. 

The movement of the steep negative wave along 
the reservoir accelerates water in the reservoir. 
Many factors affect this acceleration so that the 
approximations of the true outflow hydrograph 
depend on the factors either taken into account or 
neglected. The accelerated discharges are considered 
here as a fictitious inflow hydrograph to be taken as 
given either when only the storage continuity equa- 
tion P — Q=dS/dT is used for routing procedure or 
when both the continuity and dynamic equations 
are applied. The degree of unsteadiness of water 
movement is the criterion for the use of only the con- 
tinuity equation or of both it and the momentum 
equation for the routing. The greater the height of 
the steep negative wave in relation to the reservoir 
depths, the higher and more unsteady the inflows, 
and the greater the flow resistance, then the higher 
will be the degree of unsteadiness of flow motion, 
and the use of the two de Saint Venant equations 
with the computation involved will be more neces- 
sary. The development of expressions for the quan- 
titative determination of this criterion is a special 
problem not treated here. A criterion for the effect 
of flow resistance on the movement of steep negative 
wave is given by eq (28). When the ratio e is higher 
than a selected limit, the flow resistance effect has 
to be taken into account, regardless of the routing 
procedure. 

The fact that the steep wave is a discontinuity, 
not suitable to be considered as an unsteady flow 
gradually varied, and the fact that the movement 
and transformation of steep negative wave is not 
affected by small downstream disturbances in the 
reservoir, make the treatment of sudden water 
release different from an ordinary routing procedure. 
The determination of the outflow hydrograph can 
be divided into two procedures: (a) The effect of 
steep negative wave and of flow resistance on the 
water acceleration, combined with the inflow hydro- 
graph during the period of wave movement; (b) 
the routing procedures, based on the position in 
place and time of accelerated discharges and inflow 
discharges along the reservoir, as boundary conditions. 
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